Euler buckling instability and enhanced current blockade in suspended single-electron transistors 
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Single-electron transistors embedded in a suspended nanobeam or carbon nanotube may exhibit effects orig- 
inating from the coupling of the electronic degrees of freedom to the mechanical oscillations of the suspended 
structure. Here, we investigate theoretically the consequences of a capacitive electromechanical interaction 
when the supporting beam is brought close to the Euler buckling instability by a lateral compressive strain. 
Our central result is that the low-bias current blockade, originating from the electromechanical coupling for the 
classical resonator, is strongly enhanced near the Euler instability. We predict that the bias voltage below which 
transport is blocked increases by orders of magnitude for typical parameters. This mechanism may make the 
otherwise elusive classical current blockade experimentally observable. 

PACS numbers: 73.63.-b, 85.85.+j, 63.22.Gh 



I. INTRODUCTION 

Single-electron transistors (SET) are extremely sensitive 
devices, and are investigated as position detectors for nano- 
electromechanical systems (NEMS).'~ 3 But the reduced size 
of the mechanical resonator implies that the back-action of 
the SET can have significant effects on the mechanical degree 
of freedom, such as the generation of self-oscillations. 4 6 In 
practice, the detector and the resonator have to be investigated 
collectively as a single system. A prominent example of the 
new effects displayed by this device is the current blockade 
that appears at low bias voltage V when the SET is coupled 
capacitively to a classical oscillator. 7 The physical idea be- 
hind this phenomenon is simple: The presence of an extra 
electron on the central island of the SET induces an additional 
electrostatic force F e on the oscillator (see Fig. 1). The equi- 
librium position of the oscillator is thus shifted by a distance 
F e /k, where k is the oscillator spring constant. After such a 
displacement, the gate voltage V g seen by the SET changes by 
a quantity of the order of F e x FJek = E-^je, where e is the 
electron charge. The dimension of the conducting window in 
V s is controlled by V, since at low temperatures current can 
flow through the device only if |V g | < V (when measuring 
V g from the degeneracy point). Thus, for eV < Ee, the fluc- 
tuation of the electronic occupation of the central island suf- 
fices to bring the device out of the conducting window. The 
current is blocked for eV < E E and a mechanical bistability 
appears. 7 This phenomenon is the classical counterpart of the 
Franck-Condon blockade in molecular devices 8,9 that has re- 
cently been observed in suspended carbon nanotubes for high- 
energy vibrational modes. 10 The classical case has been theo- 
retically studied in the case of a single-level quantum dot, 11 ' 12 
as well as in the metallic case. 7 ' 13-15 

Recent experiments 1617 on suspended carbon nanotubes 
have observed a reduction of the mechanical resonance fre- 
quency of the fundamental bending mode at low bias voltages 




gate 
V„ 



FIG. 1. Sketch of the considered system: a suspended doubly- 
clamped beam forming a quantum dot electrically connected to 
source and drain electrodes held at chemical potentials p L and p R 
by the bias voltage V, respectively. The beam is capacitively coupled 
to a metallic gate kept at a voltage V g , which induces a force F c that 
attracts the beam towards the gate electrode. An additional, exter- 
nally controlled compressional force F acts on the beam and induces 
a buckling instability. 



and for V g near the degenerate region. This effect is a precur- 
sor of the mechanical instability and thus of the current block- 
ade. But the complete observation of the latter phenomenon 
is difficult since the typical value of Ee is only of a few yueV, 
thus smaller than cryogenic temperatures. In order to increase 
Ee, one can increase the electrostatic coupling between the 
oscillator and the SET since E E depends quadratically on F e . 
But another way of strengthening the effect would be to re- 
duce the spring constant k of the oscillator. The reason is that 
softer oscillators will displace more under the influence of the 
electrostatic force F e , and thus will see a larger change in the 
gate voltage when electrons tunnel in. A way of reducing k 
in a controlled manner is to operate a doubly-clamped beam 
subject to a lateral compression force F. The latter can bring 
the beam to the well-known Euler buckling instability 1 8 ' 19 (see 
Fig. 1). Under the action of the force F, the system exhibits a 
continuous transition from a flat to a buckled state, while the 
fundamental bending mode becomes softer as one approaches 
the mechanical instability (k — > 0). It is clear that this does 
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not imply a divergence of Ee, since at the transition anhar- 
monic terms will modify the simple arguments given above. 
However, a strong enhancement of E E is expected. The Eu- 
ler instability has been studied both experimentally 20-23 and 
theoretically 24-27 in micro- and nanomechanical systems. We 
have recently considered the Euler instability in NEMS for 
the case where F e is negligible with respect to the intrinsic 
electron-phonon coupling. 28 

In this paper we investigate in detail the idea of increas- 
ing the current blockade by exploiting the Euler instability, 
considering how the anharmonic terms, the temperature, and 
the non-equilibrium fluctuations modify the simplified picture 
given above. We find that near the buckling instability, the cur- 
rent blockade induced by the mechanical resonator is strongly 
enhanced, rendering this effect experimentally observable. 

The paper is structured as follows: In Sec. II the model used 
to describe the system is introduced. In Sec. Ill, a statistical 
description in terms of a Fokker-Planck equation is given. It 
is then used in the remainder of the paper to determine the 
current and the mechanical behavior of the system. In Sec. IV, 
the enhancement of Ee is obtained at mean-field level. We 
discuss the effects of thermal and charge fluctuations on the 
results in Sec. V. In Sec. VI, we investigate the consequences 
of a finite average excess charge on the quantum dot for our 
results. Sec. VII presents some estimates of the effect we are 
predicting for recently-realized experiments. We conclude in 
Sec. VIII. We relegate several technical issues to appendices. 



H. MODEL 

As a representative model for the problem outlined in 
the Introduction, we consider a quantum dot embedded in a 
doubly-clamped beam as shown in Fig. 1 . The presence of the 
metallic gate near the dot is responsible for the coupling of the 
bending modes of the beam to the charge state of the dot. The 
Hamiltonian of the system can then be written as 

H = ffyib + #SET + H c , (1) 

where // V ; D describes the oscillating modes of the nanobeam, 
//set the electronic degrees of freedom of the single-electron 
transistor, and H c the coupling between the SET and the res- 
onator. The model describes, for instance, transport through 
suspended carbon nanotubes as considered in the experiments 
of Refs. 10, 16, 17, and 29. Notice that the model also de- 
scribes an alternative setup that may be experimentally re- 
alized, namely a non-suspended quantum dot coupled to a 
beam-like gate electrode to which a compressive strain is ap- 
plied. 

Using standard methods of elasticity theory one can show 
that, close to the buckling instability, the frequency a> of the 
fundamental bending mode of the nanobeam vanishes while 
those of the higher modes remain finite. 19 This allows one to 
retain only the fundamental mode parametrized by the dis- 
placement X of the center of the beam. As detailed in Ap- 
pendix A, the Hamiltonian representing the oscillations of the 



nanobeam thus takes the Landau-Ginzburg form 21 ' 24 " 26 ' 28 

P 2 into 2 9 a a 

# vib = — + X 2 + -X\ (2) 

2m 2 4 ; 

where P is the momentum conjugate to X. For a doubly- 
clamped uniform nanobeam of length L, linear mass density 
cr, and bending rigidity k, one can show 25 ' 26 that close to the 
instability the effective mass of the beam is m — 3crL/8. The 
fundamental bending mode frequency reads 



w = w yi-— , (3) 

where F is the compression force, F c = K(2n/L) 2 the critical 
force at which buckling occurs, and ojq = v 'k I 'cr(2n / L) 2 . The 
positive parameter a - F c L(n/2L) 4 ensures the stability of 
the system for F > F c . 30 For F < F c (w 2 > 0), X = is 
the only stable solution, and the beam remains straight. For 
F > F c , it buckles into one of the two metastable states at X = 
+ yj-mcj 2 ja. Notice that in writing Eq. (2), we assumed that 
the nanobeam cannot rotate around its axis due to clamping at 
its two ends. 

Electronic transport is accounted for by the SET Hamilto- 
nian consisting of three parts, 

HSET = H Aot + //leads + Stun. (4) 

where ff<iot describes the quantum dot, //i ea ds the left (L) and 
right (R) leads, and H tun the tunneling between leads and dot. 
Explicitly, 

/ — \ U 
//dot = (ed -eV g ) « d + — «d(«d - 1), (5) 

with n& = (fid, and <fi (d) creates (annihilates) an electron on 
the dot, V g = CgVg/Cs, with C g and Cj; the gate and total 
capacitances of the SET, respectively. The intra-dot Coulomb 
repulsion is denoted by U. In the following, we set e& = 0, 
measuring V g from the degeneracy point. The left and right 
leads are assumed to be Fermi liquids at temperature T with 
chemical potentials /iL and ^ur (measured from ed), respec- 
tively. A (symmetric) bias voltage V is applied to the junction 
such that jUl = _ A(r = eV/2. The lead Hamiltonian reads 

Hi eais = J](e k -fi a )cl a c ka , (6) 

ka 

with q„ the annihilation operator for a spinless electron of 
momentum k in lead a = L, R. 31 Finally, tunneling is ac- 
counted for by the Hamiltonian 

//tun = Yl {taC la d + h - C - } ' (?) 

ka 

with t a the tunneling amplitude between the quantum dot and 
lead a. 

Two different kinds of couplings exist between the elec- 
tronic occupation of the dot «d and the vibrational degrees of 
freedom: (i) an intrinsic one that originates from the varia- 
tion of the electronic energy due to the elastic deformation 
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of the beam, and (ii) an electrostatic one, induced by the 
capacitive coupling to the gate electrode of the SET. 33-36 By 
symmetry, the former is quadratic in the amplitude X and its 
effect on the Euler instability has been considered in Ref. 28. 
The latter is linear in X and here we are interested in the case 
where the second coupling dominates over the first one. Their 
relative intensity is controlled by the distance h between the 
gate electrode and the beam, since the intrinsic coupling does 
not depend on h, while the electrostatic force depends log- 
arithmically on h. 33 Assuming that the beam is sufficiently 
close to the gate electrode such that the capacitive coupling 
dominates, 37 we can write 

ff c = F e Xn A , (8) 

where — F e is the force exerted on the tube when one excess 
electron occupies the quantum dot (see Fig. 1). The model as- 
sumes that the gate voltage is such that only charge states with 
«d = and 1 are accessible. For larger gate voltages overcom- 
ing the charging energy of the quantum dot, the charge on the 
dot will instead fluctuates between N and N + 1 . This induces 
an additional constant force bending the tube further. The ef- 
fect of such a force on the classical current blockade will be 
discussed in Sec. VI. Notice that in the case the suspended 
structure is the gate capacitance coupled to a static quantum 
dot, the intrinsic coupling is not present, and only the capaci- 
tive electromechanical coupling has to be taken into account. 

IIL FOKKER-PLANCK DESCRIPTION 

We are interested in describing the vicinity of the instability 
where the relevant resonator frequency vanishes [see Eq. (3)]. 
The mechanical degree of freedom can then be treated classi- 
cally since for any reasonable temperature, ftco <sc k B T. The 
softening of the mechanical mode implies also a natural sep- 
aration of timescales between the slow mechanical mode and 
the fast electronic degrees of freedom, controlled by the typ- 
ical tunneling rate F. As detailed in Appendix B, it is thus 
convenient to eliminate the fast modes and obtain a Fokker- 
Planck equation for the probability distribution P(X, P, t) of 
the slow mode: 7 - 1 u4 - 28 - 38 

d,<P = - -d x P - F &s (X)d P P + ?7(X) + ,??e dp(PP) 
m m 

J^+^tW (9) 

The effective force F a s(X) = -dxH^b + F c .i(X) acting on 
the mechanical degree of freedom consists of two parts: a 
force arising from the Hamiltonian (2) of the nanobeam, 
-dxH v n, = -mcj 2 X - aX 3 , and a current-induced conserva- 
tive force F c -i(X) = -F e no(X), proportional to the occupa- 
tion of the dot averaged over a time long with respect to r~ , 
but short with respect to the period of the mechanical mo- 
tion, no(X) = («d)x- In Eq. (9), the diffusion constant D(X) 
accounts for the fluctuations of the force associated with the 
coupling Hamiltonian (8) originating from the stochastic na- 
ture of the charge-transfer processes. Finally, retardation ef- 



fects cause dissipation of the mechanical energy, with damp- 
ing coefficient rj(X). 

To account for the quality factor Q - muolr\ & of the 
nanobeam mode, the mechanical degree of freedom is cou- 
pled to an additional environment at equilibrium (such as, e.g., 
a generic phonon bath within the Caldeira-Leggett model 39 ), 
implying dissipation and fluctuations controlled by an extrin- 
sic damping constant entering Eq. (9). This extrinsic damp- 
ing comes from several mechanisms coupling the mechani- 
cal mode to other degrees of freedom: localized defects at 
the surface of the sample (thought to be the main source of 
dissipation in semiconductor resonators 40 and which can be 
modeled as two-level systems 41,42 ), clamping losses, thermo- 
elastic losses, ohmic losses due to the gate electrode (which 
have been predicted to be the dominant source of extrinsic dis- 
sipation for graphene-based resonators in Ref. 43), etc. Due 
to the wide variety of these possible sources of extrinsic dis- 
sipation, we here assume for simplicity that they can all be 
lump into the generic (Ohmic, memory-free 39 ) damping con- 
stant r] e . Notice also that the phonon temperature of the bath is 
typically lower, but of the same order as the electronic temper- 
ature T. 16 For simplicity, we assumed in writing Eq. (9) that 
both temperatures coincide, as we do not expect a qualitative 
change of our results due to this assumption. 

The explicit form of the coefficients entering into the 
Fokker-Planck equation (9) depends on the transport regime 
one considers (sequential or resonant transport), as well as on 
the nature of the quantum dot (metallic or single-level quan- 
tum dot). In this paper we consider the case of a single level in 
the sequential tunneling regime, but a similar analysis can be 
carried out for the metallic (e.g., along the lines of Refs. 7, 14, 
28) and the resonant transport regime (cf. Refs. 11, 12, 38). 
To be specific, we assume that HT - 2 fl =L,R KT a <sz k^T, with 
Y a — 27r\t a \ 2 v/fi and v the density of states at the Fermi level 
of the leads. We also assume the intra-dot Coulomb repulsion 
U — ♦ °o such that double occupancy of the dot is forbidden. 
In this transport regime, the position-dependent rates for tun- 
neling into and out of the dot read 44 

vn (F e X~eV g -fi a \ 

r 0l (x) = V r fl / F g , (io) 

^ [ (FvX-eVe-UaX 

r w (X)= V r a i-/ F , (ii) 

respectively, where fp(z) = (e z + 1) _1 is the Fermi function. 
Thus, the average occupation of the dot for a given mode am- 
plitude X is 

™ r Q1 (X) 

no(X) = , (12) 

and, as shown in Appendix B, we have D(X) = 2F 2 «o(X)[l - 
«oPO]/r and rj{X) = -F e dxno{X)/T for the current-induced 
diffusion and damping terms in Eq. (9). 14 ' 45 The average cur- 
rent / through the device can be obtained from the stationary 
solution of the Fokker-Planck equation (9), d/P^ = 0, by av- 
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eraging the position-dependent current 

(F e X - eV g - eV/2 



I(X) 



^r L r R 



h 



-/f 



F e X- 



eV/2 



k R T 



with the phase-space distribution, 

! = jj dXdPP st (X,P)I(X). 











,1 










20 - 


s *t — 

V / 
/ • ( 




(13) 




-H 1 










U. 

-1 


-0.5 
1 




1 1 1 =f 

\ 0.5 

\. s 


1 > 

l 6 


(14) 






-20 > 







Before we proceed, it is convenient to introduce reduced 
variables in terms of the relevant energy scale of the problem 
£g = Fg/moiS, the polaronic shift £ = Fe/maJi., and the vi- 
brational frequency for vanishing compression force ojq [see 
Eq. (3)]. By denoting x — X/{, p — Pjrrmot, t = a>^t, the 
Fokker-Planck equation (9) becomes 



8 T P = - pd x r - f ett (x)d p <P + (y(x) + y e )8 p (pP) 

with the scaled effective force given by 

fes(x) — 6x — ax 1 — uq(x), 



(15) 



(16) 



the reduced force <5 = F/F c — 1, and the anharmonicity pa- 
rameter a = a£ 4 /E^. We further introduced a scaled current- 
induced diffusion constant 



2u)q 

d(x) = -—n (x) [1 - n (x)] , 



and damping coefficient 



y(x) = -^r3 x «o(x). 



(17) 



(18) 



In Eq. (15), y e = rj e /mu>o = Q~ l , where Q is the quality factor 
of the mechanical resonator and T = k^T In these scaled 
units, the electromechanical coupling appears only in the co- 
efficient of the quartic term, a = aF\ / Xmafy 3 . It is impor- 
tant to notice that for actual experiments on suspended carbon 
nanotubes, 16 ' 17 ' 29 5 « 1 as we will discuss more extensively 
in Sec. VII. 



IV. MEAN-FIELD APPROACH: ENHANCEMENT OF THE 
CURRENT BLOCKADE 

We begin our analysis by assuming wq/F — > 0. 46 Note 
that the diffusion and dissipation coefficients d(x) and y(x) 
in Eq. (15) are proportional to wo/E [cf. Eqs. (17) and (18)], 
so that this implies to neglect current-induced fluctuations. In 
this limit, the stationary solution for f is given by a Boltz- 
mann distribution at temperature f , 



^st(x, p) = N exp 



p 2 /2 + v eS (x) 



(19) 



FIG. 2. Example of a solution x{6) of the equation for dynamic 
equilibrium (21) for rc - (dashed line), n = 1 (solid line), and for 
no = 1 and a = (dotted line). In the figure, a = 10~ 3 . 



with N a normalization constant. In order to obtain trans- 
parent analytical results, we also assume zero temperature (in 
fact, « ^ B r « Eg), such that the stationary probability 
distribution (19) becomes r P st (x,p) = 6(p)6(x- x m ). Here, x m 
is the global minimum of the effective potential 



v e «(x) = - J dx'f eS (x') 



(20) 



corresponding to the effective force (16) and can be deter- 
mined from the dynamical equilibrium equation 



f e s(x) = 0, 



d/eff(x) 

dx 



< 0. 



(21) 



Notice that the latter equation can have more than one so- 
lution, such that the system is multi-stable. In this zero- 
temperature limit, the current can then easily be obtained from 
Eq. (14). Doing so as a function of the gate and bias voltages, 
we can determine the Coulomb diamond for a given compres- 
sion force 5. At zero temperature, one finds that there always 
exists a region at low bias voltage where the current is sup- 
pressed. 

To characterize this classical current blockade, we define 
A„, the minimal value of bias voltage, for which a finite cur- 
rent flows through the device at zero temperature. It is useful 
to first derive a simple estimate of the maximally obtainable 
A„. To do so, we solve the dynamic equilibrium equation (21) 
for «o = and «o = F corresponding to empty and occupied 
central island, respectively. For «o - one has the solutions 
x — for any 6 and x = + y/6/a for 5 > of the pristine Euler 
instability (see Fig. 2, dashed line). For «o = 1, the solutions 
can easily be sketched for a <k 1 as an interpolation of the 
a — solutions (dotted line in Fig. 2) and the solutions for 
«o = 0. The exact result is shown as a solid line in Fig. 2. 
We are interested in the maximum shift in x that the system 
undergoes in response to a fluctuation of no by one unit, Ax. It 
is apparent from the figure that this happens for 6 — where 
Ax = The corresponding change in the effective po- 

tential (20) is Au ef f ~ 1 / This provides an estimate of 
the maximal energy gap generated by the electromechanical 



coupling, and thus a good estimate of A„. Notice that the sim- 
ple argument above is not specific to the transport model we 
are considering here, as the specific form of «o in the con- 
ducting region does not enter our argument. Thus, we expect 
that our estimate of a maximal gap A„ ~ 1 / "\f& remains valid 
for metallic quantum dots, as well as in the resonant transport 
regime. 

We now turn to the complete solution of Eq. (21). For sim- 
plicity, we assume symmetric coupling to the leads (Tl = 
Tr = F/2), such that the average occupation of the dot at fixed 
x [entering into the effective force (16)] is obtained from the 
zero-temperature limit of Eqs. (10) and (12) and given by 



n (x) = 



1 



(22) 



where v = eV/Zig (assumed positive for definiteness), v s = 
eVg/fig, and Q(z) is the Heaviside step function. In the most 
general case, we solve Eq. (21) numerically. However, trans- 
parent analytical expressions can be obtained in the limits 
\S\ » \f& and \5\ <K if&. In particular, we can obtain explicit 
expressions for the value of v beyond which the current begins 
to flow, i.e., the gap A„. To first order in the small parameter 
yf&/\6\ (far from the instability) and \6\/ if& (in the vicinity of 
the instability), we find that 



1 

~26' 
1 

46' 

1/2-1 1 3 

lf& 24/3 



- 1 < 6 <sc - 



6 » if&, 



(23) 



Far from the mechanical instability, the gap is simply given 
by the result of a harmonic theory (see Appendix C) where 
A„ = 1 /2u"(x), with v"(x) the curvature of the bare potential 
v(x) [i.e., the effective potential without the contribution from 
no(x)] at its global minimum x. Far below the buckling insta- 
bility (-1 < 6 <sc -lf&), x = and v"(x) = -6, such that 
A„ = — 1 125. Far above the mechanical instability (5 » V&), 
x - - y[6~J& and v"(x) = 26, such that A„ = 1/45. As one 
approaches the buckling instability from below or above, the 
apparent divergences in the first two lines of Eq. (23) are cut- 
off by the cubic term in x in the effective force (16), and for 
\6\ <K y&, the maximal gap A„ ~ 1/ y& is reached. 

The analytical results of Eq. (23) are compared to a nu- 
merical calculation of the gap in Fig. 3(a) for a = 10~ 3 and 
a = 1CT 6 (red dots and blue squares in the figure, respec- 
tively). It is evident from the figure that there is a dramatic 
increase of the gap close to the instability. Furthermore, the 
smaller a, i.e., the smaller the electromechanical coupling, 
the larger is the increase of the gap at the instability relative 
to its value for vanishing compression force [see the inset in 
Fig. 3(a)]. However, of course, the maximal value of the gap 
in absolute terms increases with the strength of the electrome- 
chanical coupling as F A J 2 . It would thus be of great experi- 
mental interest to exploit the Euler instability to obtain a clear 
signature of the classical current blockade in transport exper- 
iments on suspended quantum dots. 
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FIG. 3. (Color online) (a) Gap A„ and (b) gate voltage v g (defined as 
the bias and gate voltages in reduced units at the apex of the Coulomb 
diamond, respectively) as a function of the scaled compression force 
6 = F/F c - 1. The red circles and blue squares are numerical re- 
sults for a = 10~ 3 and 10~ 6 , respectively, which are compared to 
the asymptotic behaviors (23) and (24) for forces below (solid line), 
above (dashed line), and in the vicinity (dotted line) of the critical 
force F c . Inset: Gap A„ ~ 1/ from Eq. (23) as a function of & at 
the mechanical instability (S = 0). 



The gaps of Eq. (23) are obtained for values of the gate 
voltage approximately given by 



1 

25" 

- 1 - F 

A6 V ff' 

_L (3 + 



-1 < 5 <sc -if&, 



6 » -ya, 



\6\ <sc if&, 



(24) 



which are shown in Fig. 3(b) and compared to a numerical 
calculation. Equations (23) and (24) define the apexes of 
the Coulomb diamonds which are shown in Fig. 4. The ef- 
fect of the compression force is thus to continuously displace 
the Coulomb diamond in the v-u g plane towards negative gate 
voltages [see also Fig. 3(b)], and to open a gap which is maxi- 
mal close to the Euler instability at 6 = [see Fig. 4(d)]. Note 
that the shift in gate voltage is strongly asymmetric about the 
Euler instability. While the shifts are only small below the 
Euler instability [see Figs. 4(a)-(c) and Fig. 3(b)], the shifts 
in gate voltage are orders of magnitudes larger on the buckled 
side of the Euler instability [see Figs. 4(e)-(g) and Fig. 3(b)]. 
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FIG. 4. (Color online) Mean-field current / at zero temperature and for symmetric coupling to the leads (T L = T R = T/2), as a function of bias 
v and gate voltage v g (measured in units of the elastic energy The (scaled) compression force 6 increases from (a) to (g). Notice that the 
scale of the u g -axis is different in (e), (f), and (g), and in (a)-(d). The red dashed lines indicate the position of the Coulomb diamond in the 
absence of electromechanical coupling (F c = 0). In the figure, a = 10~ 6 , and dark blue and white regions correspond to / = and / = eT/4, 
respectively. 



In fact, it may be that these shifts would be the most easily de- 
tected consequence of the Euler buckling instability in NEMS. 
In Fig. 4, the bias and gate voltages are measured in units 
of the elastic energy which is of the order of a few ytzeV 
for typical experiments on suspended carbon nanotubes (see 
Sec. VII). The smallness of this energy scale explains why the 
scaled numerical values of the shifts become so large on the 
buckled side of the Euler instability. 

It is also interesting to comment on the shape of the 
Coulomb blockade diamond. In Ref. 28, we showed for the 
case of intrinsic electron-phonon coupling (quadratic in x) and 
for a metallic quantum dot that the Euler buckling instability 
leads to nonlinear deformations of the Coulomb diamonds, 
a phenomenon that we have named "tricritical current block- 
ade". In contrast, our present results show that for a capacitive 
electromechanical coupling (linear in x) and for a single-level 
quantum dot, the shape of the Coulomb diamond remains un- 
changed. The conventional triangular shape in the v-v g plane 
is delineated by straight lines with u ~ +2u g for any value 
of the compressive strain. As we have checked, 47 the differ- 
ence between the present results and those of Ref. 28 is due 
to the difference in the transport models considered (metal- 
lic vs. single-level quantum dot), and not to the type of elec- 
tromechanical coupling (intrinsic vs. extrinsic). Specifically, 
we find that the difference is due to the fact that in the single- 
level case, the average occupation of the dot abruptly jumps 
as a function of gate voltage, while in the metallic case, this 
occupation gradually changes due to the continuous density 
of states of the dot. Notice also that for intrinsic electrome- 
chanical coupling, the Coulomb diamond is not influenced 



by the latter for compression forces below the critical force 
(5 < 0), in contrast to the present case. This is due to the fact 
that, in the flat state, the quadratic electromechanical coupling 
merely represents a renormalization of the fundamental bend- 
ing mode frequency, and does not lead to current blockade. 

It is instructive to make the analogies with standard re- 
sults of Landau mean field theory for continuous phase 
transitions 48 explicit. According to Eq. (21) governing the 
dynamical equilibrium, we can make the following identifica- 
tions: x corresponds to the order parameter in Landau theory, 
6 to the reduced temperature, and a to the coefficient of the 
quartic term in the Landau free energy. Finally, «o plays a 
role similar to a symmetry-breaking (magnetic) field. In the 
present case, this field is in general dependent on x, which 
has no correspondence in Landau theory. Nevertheless, the 
analogy between no and a magnetic field is helpful since some 
of our results can be understood by comparing the situations 
with zero (no - 0) and one (no = 1) electrons on the dot, 
as illustrated by the above estimate for the maximal A„ (see 
Fig. 2). 

With these correspondences, we can now establish analo- 
gies between some of our results and standard results of Lan- 
dau theory. To start with, the dependence of the displacement 
x ~ ±6 1 ^ 2 in the buckled state is analogous to the result of 
Landau theory that the order parameter exponent is J3 = 1 /2. 
Given that A„ depends linearly on x, we can also interpret the 
relations in Eq. (23) in terms of Landau theory. Let us start 
with the case of small 6, in the immediate vicinity of the in- 
stability. In this case, we find that A„ ~ a -1 ^ 3 . The exponent 
of a corresponds to the critical exponent 5 — 3 of Landau the- 
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ory governing the dependence of the order parameter on the 
symmetry-breaking field at the critical temperature. Further 
from the instability, we have A„ ~ 1/|<S|. This relation is re- 
lated to the familiar Curie law for the order parameter (or the 
susceptibility) as function of temperature in an external field, 
with mean-field critical exponent y = 1. 

Let us finally emphasize that the analogy with Landau the- 
ory is restricted to mean-field level, since contrarily to critical 
phenomena where an infinite number of modes is present, the 
system we describe is constituted by a single mode. Moreover, 
beyond mean-field theory, fluctuations in Landau theory are 
purely thermal, while in the present context, non-equilibrium 
fluctuations play an essential role. It is the effects of these 
fluctuations which we turn to in the next section. 
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V. THERMAL AND CURRENT-INDUCED 
FLUCTUATIONS 

We now go beyond the mean-field results of the previous 
section by taking into account the effects of the thermal as well 
as the current-induced fluctuations. It is physically clear that 
these fluctuations will lead to a smoothening of the current 
blockade at low bias voltages, as the system can explore more 
conducting states in phase space. 



A. Temperature effects 

We first neglect the current-induced fluctuations and focus 
on thermal fluctuations only. As discussed in Sec. IV, this 
becomes asymptotically exact in the extreme adiabatic limit 
of u>o/r — > 0, where the terms y(x) and d(x) can be dropped 
from Eq. (15). The stationary solution for f is then given by 
the Boltzmann distribution (19) with the effective potential 



5x 2 ax 4 T I ( x - v g + v/2 

v eS (x) = - — + — + x + - In I f F 



T I lx-v g -vl2 



(25) 



The current can now be easily calculated by numerical inte- 
gration of Eq. (14) with Eq. (19). The result is shown in Fig. 5 
as a function of the bias voltage, for gate voltages correspond- 
ing to the apex of the modified zero-temperature Coulomb 
diamond [cf. Eq. (24)]. Once plotted as a function of v/A v , 
one finds that the current behavior is similar at the transition 
(Fig. 5) and far from the transition (inset of Fig. 5). In both 
cases, the low-bias blockade of the current becomes less pro- 
nounced as temperature increases, and vanishes completely 
for temperatures of the order of the gap A„. As shown in Ap- 
pendix D [cf. Eq. (Dl)], the current has a Fermi-function-like 
behavior as a function of the bias voltage for temperatures 
much smaller than the energy gap A„ (see dashed and dashed- 
dotted lines in Fig. 5). It is thus exponentially suppressed for 
bias voltage below the gap. At larger temperatures, Eq. (D4) 
shows that the current is linear in the bias voltage (see dotted 
line in Fig. 5). 



FIG. 5. (Color online) Current / at the apex of the Coulomb diamond 
as a function of bias v scaled by the energy gap A„ for various values 
of r/A„, and for compression forces in the vicinity of the buckling 
instability (\6\ <k yf&). In the figure, only the temperature-induced 
fluctuations are considered. Inset: Same as the main figure for com- 
pression forces far from the buckling instability (|<5| » ~\f&). 



Our numerical and analytical results (cf. Appendix D) thus 
confirm that tuning the system near the buckling instability 
where A„ dramatically increases allows one to enlarge the tem- 
perature region over which the current blockade is observable. 



B. Non-equilibrium dynamics close to the mechanical 
instability 

We now consider the non-equilibrium Langevin dynamics 
of the nanobeam by solving the full Fokker-Planck equation 
(15). This is done by discretization of the Fokker-Planck 
equation and solution of the resulting linear system. We fo- 
cus on the transition region (S = 0) and calculate the current 
for y g at the apex of the Coulomb diamond [see Eq. (24) and 
Fig. 3(b)] and temperature lower than the gap f/A„ = 0.1. Be- 
fore we present our results, we notice that for (cjq/T, y e ) «: 1, 
we can show that the stationary distribution of the Fokker- 
Planck equation approximately only depends on the ratio 
a result we have also checked numerically (see Appendix E for 
details). The reason for this behavior is that, for (o>q/T, %,) <K 
1, the stationary distribution is almost a function of the (re- 
duced) energy E = p 1 12 + v e s(x) only. 

Numerical results for the current are shown in Fig. 6 for var- 
ious ratios of the inverse quality factor Q l = y e as quantified 
by the damping coefficient y e and the adiabaticity parameter 
o>o/r. Our principal observation is that the current blockade 
becomes sharper for low-Q resonators. 

One can qualitatively understand the behavior of the current 
in Fig. 6 by defining an effective temperature of the system 



Teff - 



(d)/2 + y e f 

(j) + 7e 



(26) 



in close analogy with the fluctuation-dissipation theorem. 49 
In Eq. (26), (d) and (y) are the averages over the phase- 
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FIG. 6. (Color online) Current / at the apex of the Coulomb dia- 
mond for 6 = as a function of v/A D for a = 10~ 6 and f/A„ = 0.1. 
Solid lines: non-equilibrium Langevin dynamics for = 1, 10" 1 , 
10~ 2 , 10~ 3 , and 0, from the highest to the lowest curve at large 
bias. Dashed-dotted line: fully adiabatic limit = oo), i.e., cur- 
rent only including thermal fluctuations (cf. dashed-dotted curve in 
Fig. 5). In our numerical calculations, we used a>o/T = 10~ 2 . 



space probability distribution of the current-induced fluctu- 
ations and dissipation [cf. Eqs. (17) and (18)], respectively. 
Notice that the strength of these two quantities is controlled 
by the adiabaticity parameter wq/F. 

As one can see from Fig. 6, for u < A„, the current is almost 
insensitive to the quality factor, and is the same as without 
current-induced fluctuations (see dashed-dotted line in Figs. 5 
and 6). Using Eqs. (17) and (18), we have 



d(x) = 2y(x)T + 



wo 
2Y 



/f 



v/2 



/f 



■ Va + V/2 



T I \ T 

(27) 

for symmetric coupling to the leads. However, for v < A„, 
positions x for which the current I{x) of Eq. (13) is sup- 
pressed are most stable (see dashed line in Fig. 7), such 
that the current-induced diffusion and damping constants ap- 
proximately satisfy a "local" fluctuation-dissipation theorem 
for all relevant positions x that are significantly populated, 
d(x) ^ 2Ty(x). We thus have (d) ^ 2T(y), and according to 
the definition (26), we have T e ff - T . Hence, for bias voltages 
lower than the energy gap A„, the current-induced fluctuations 
behave as the thermal ones, essentially keeping the mechani- 
cal system at equilibrium. 

On the contrary, for v > A„, positions x for which the sys- 
tem is conducting are the most stable ones (see dotted line in 
Fig. 7), and one has (d) - <wo/2F, while (y) is exponentially 
small. The mechanical system is then subject to strong non- 
equilibrium fluctuations. For (y) <k y e , we thus have from 
Eq. (26) r e ff — f + . This estimate of the effective tem- 
perature shows that the system becomes "hotter" as the ratio 
decreases. 50 Hence, the system can explore more states 
in phase space for which I(x) is suppressed, and, in turn, the 
current decreases for decreasing f° r u > A„ (see Fig. 6). 
The latter argument breaks down when %, ^ il)- In tnat case, 
we can estimate the effective temperature self-consistently, by 
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FIG. 7. (a) Current J as a function of x [Eq. (13)] and (b) effec- 
tive potential D e ffM [Eq. (25)] for bias voltages below (dashed line), 
above (dotted line), and at (solid line) the energy gap A„. The param- 
eters are the same as in Fig. 6, i.e., 6 = 0, a = 10~ 6 , 77 A„ = 0.1, and 
v g = -3/4^1. 



assuming that the phase-space distribution is a Boltzmann dis- 
tribution at the temperature T e s, Approximating the effective 
potential by its zero temperature expression, and averaging 
d(x) and y{x) over the effective Boltzmann distribution, we 
find for u » A„ 



nA D 
128A 



exp 



Am 1 



(28) 



with A = 9(1 - 2- 1/3 )/2 11/3 . We thus have f eff /A„ ~ 
(v/A„) 2 / In (v/A v ) » T/A„, which explains why for y e = 
the current is more suppressed than for finite y e . It is also in- 
teresting to note that this estimate of the effective temperature 
is much larger than for a metallic quantum dot, where T e g ~ v 
(Ref. 13). The reason for this difference is that, in the metallic 
case, the fluctuation and the dissipation are of the same order 
inside the bias window, while in the single-level case, the av- 
erage dissipation is exponentially suppressed as y(x) only has 
a significant contribution for positions x corresponding to the 
borders of the Coulomb diamond [see Eq. (18)]. 

Our results show that a low quality factor is more suitable 
for the observation of the current blockade in classical res- 
onators. It is interesting to note that this conclusion is also 
valid in the quantum case, 8 ' 9 where the Franck-Condon block- 
ade is more pronounced for fast equilibration of the vibron 
mode. Due to the scaling of our results for the classical cur- 
rent blockade with the parameter (see Fig. 6), we also 
conclude that it is advantageous for the observation of this 
phenomenon to have a resonator which is slow compared to 
the tunneling dynamics, i.e., a>o «: T. 
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VI. EFFECT OF A FINITE EXCESS CHARGE ON THE 
QUANTUM DOT 

Within the transport model of a single resonant electronic 
level with infinite charging energy that we have used so far, 
the number of electrons on the dot can only vary between or 
1 (see Sec. II). More generally, the range of gate voltages can 
exceeds the charging energy and the average number of excess 
electrons N on the dot can be much larger than 1 . Due to these 
excess electrons, an additional force —Fu further bends the 
nanotube, and hence increases its vibrational frequency. 16 ' 17 
We can thus expect that the bias voltage below which the cur- 
rent is blocked will decrease when N increases. 

In order to investigate the effect of a non-vanishing average 
excess charge on the quantum dot, we assume that the gate 
voltage is such that there is either N or N + 1 electrons on the 
dot. We measure the fluctuation of the dot occupation «d with 
respect to N, and incorporate the resulting additional force in 
Eq. (16) by writing f & g(x) = Sx - ax 3 - no(x) - f N , where 
/n = FnIF s . We neglect thermal and current-induced fluc- 
tuations, and work within a mean-field approximation at zero 
temperature, such that uq{x) is given by Eq. (22). For finite 
/a?, the bare potential v(x) (i.e., without the current-induced 
contribution) can be approximated by a harmonic potential 
close to its global minimum x, such that the bias voltage below 
which the current is blocked is given by A„ = l/2v"(x) (see 
Appendix C). The energy gap (resulting from the most stable 
solution of Sx - ax 3 = fy) is plotted in Fig. 8(a), and the 
gate voltage at the apex of the Coulomb diamond in Fig. 8(b). 
As anticipated, the increase of the energy gap close to the 
mechanical instability is reduced as fy increases. Moreover, 
the displacement of the Coulomb diamond in u g is less pro- 
nounced for large fy. 

Far from (|<5| » ^Jaf N ) and in the vicinity of (\6\ <s ^Jaf^) 
the Euler instability, we analytically find for the gap 



1 

~26' 
1 

45' 

1 



A„ = 



1 - 




to first order in the small parameter ^af N /\6\ (far from the in- 



stability) and |<5|/, 



(in the vicinity of the instability). Far 



below and above the instability, the gap follows the same be- 
havior as for — [see Eq. (23)]. This is due to the fact that 
for large \S\ » /jy, the stable position of the beam is similar to 
the one for fy = 0. In the vicinity of the instability, the gap 
is reduced as fa increases as 1 jf^ 3 (see the inset in Fig. 8). 
The reduction of the maximal gap close to the instability is 
a direct consequence of the smoothening of the mechanical 
transition between flat and buckled states due to the presence 
of the symmetry-breaking force f N , similar to the behavior 
of the order parameter at a second-order phase transition in a 




FIG. 8. (Color online) (a) Gap A„ and (b) gate voltage v g at the apex 
of the Coulomb diamond for a = 10~ 6 as a function of the scaled 
compression force S for increasing values of f N [from top to bottom 
and bottom to top at S = in (a) and (b), respectively]. The blue 
squares ( f N = 0) corresponds to the numerical results of Fig. 3, while 
the solid lines result from the harmonic approximation (see text). 
Inset: Gap A„ ~ l/f^ 3 from Eq. (29) as a function of at the 
mechanical instability (S = 0). 



symmetry-breaking field. 48 

We can estimate the force above which the increase of 
the gap at the instability completely vanishes by equating in 
Eq. (29) the gap at, say, S = — 1 and 5 = 0. We obtain that 



(29) the increase of the gap should vanish once fy 

Since for large N, F^ F e N/2, this means that if the average 
charge on the dot N > 2/ V3 3 a, the increase of the gap at the 
instability completely disappears. Since a is typically small 
(see Sec. VII), we expect that a significant increase of the cur- 
rent blockade at the mechanical instability persists for a wide 
range of gate voltages. 



VII. EXPERIMENTAL REALIZATION 

The electromechanical coupling (8) is typically weak in 
experiment. For this reason, only a precursor of the classi- 
cal current blockade has been seen in two recent experiments 
on suspended carbon nanotube quantum dots, 1617 but the full 
current blockade has not yet been observed. Indeed, we can 
obtain an estimate for the frequency shift of the fundamen- 
tal bending mode, induced by the electromechanical coupling, 
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from the effective potential associated with F e ff(X). The shift 
arises from the position dependence of iiq(X). Expanding the 
current-induced force for weak electromechanical coupling, 
we find 



F c .i(X) - -F e « (0) - Fl 



dm 



deV a 



X, 



(30) 



F c =0 



i.e., the current-induced force generates a term in the effec- 
tive potential which is quadratic in X. Far from the current- 
induced instability, this term gives a small renormalization of 
the resonance frequency, Awo/<wo = (-Eg/2) dno/deV g \ F _., 
from which we can extract a reliable estimate of the energy 
scale of the current blockade, 



2C g Awo 



dn 



deV e 



(31) 



From the experiments of Refs. 16 and 17, we extract a value 
of Awo/wo of a few percents for V g at the degeneracy point. 
The derivative of no with respect to eV g can be estimated as 
the inverse of the width of the conductance peak in the V-V g 
plane, divided by Cg/Cj;. This last quantity can, in turn, be 
estimated from the slope of the Coulomb diamonds. Collect- 
ing these ingredients we find that for the suspended carbon 
nanotubes of Ref. 16, £ E - 3-5 yueV which corresponds to 
a =! 10~ 1() , while for those or Ref. 17, we get £ E 20 fieV 
and a - 10~ 8 (see Ref. 30). 

We now use these numbers to estimate the possible en- 
hancement of the current blockade near the Euler buckling 
instability. Based on Eq. (23), these parameters yield a possi- 
ble increase of the mechanically-induced gap by three orders 
of magnitude, leading to a maximal A„ (converted into a di- 
mensionful quantity using the energy scale EH) of the order 
of 3 to 5 meV. Such large gaps would be much more easily 
observable in experiment. The implementation of such a de- 
vice could be performed by the method routinely employed to 
control break junctions through a force pushing the substrate 
of the device. 

We also emphasize here again that it is preferable to operate 
the system near zero excess charge on the quantum dot, where 
there are only a few electrons on the nanotube such that the 
increase of the energy gap close to the Euler instability is not 
smeared out by the additional force exerted on the nanotube 
(see Sec. VI). However, for the parameters of Refs. 16 and 
17, we estimate that the enhancement of the current blockade 
remains very substantial for any realistic value of the excess 
charge. 

When the tunneling-induced width F becomes larger than 
or of the order of temperature, co-tunneling effects tend to 
smear the current blockade, 9 ' 12 as direct electronic transitions 
between left and right leads take place. However, close to the 
buckling instability the gap may remain larger than tempera- 
ture so that co-tunneling corrections should be suppressed in 
the immediate vicinity of the instability. 

A last comment is in order on the required precision in the 
control of the lateral compression force F. As mentioned 
above, the increase of A„ is larger for smaller a. But at the 



same time, the increase is limited to a small force range. The 
increase of the gap near the transition goes as 1/\F - F c \. 
Thus when a is very small, a stringent requirement will be 
the precision in F that we denote by 5F. In this case the max- 
imal gap will be of the order of ^ ln ^^f 1 j, assuming that 

6F/F C » if&. This result can easily be checked by convo- 
lution of the gap (23) with a Lorentzian of width 5F . This 
implies that if one is able to control the force with a precision 
sufficient to see the buckling instability (6F/F C <K 1), there 
remains a large enhancement of the gap. 



vm. CONCLUSIONS 

In this work, we have investigated the consequences of a 
capacitive electromechanical coupling in a suspended single- 
electron transistor when the supporting beam is brought close 
to the Euler buckling instability by a lateral compression 
force. Our main result is that the low-bias current blockade 
originating from the coupling between the electronic degrees 
of freedom and the classical resonator can be enhanced by 
several orders of magnitude in the vicinity of the instability. 
We show that both the mechanical as well as the electronic 
properties of this regime can be described in an asymptoti- 
cally exact manner based on a Langevin equation. These re- 
sults are a direct consequence of the continuous nature of the 
Euler buckling instability and the associated "critical slow- 
ing down" of the fundamental bending mode of the beam at 
the instability. In fact, more generally our results frequently 
have close and instructive analogies with the mean-field the- 
ory of second-order phase transitions. We focused on the 
sequential-tunneling transport regime of single-level quantum 
dots, but many of our qualitative results should remain valid 
also in the metallic case as well as for the resonant trans- 
port regime. 47 In fact, our basic approach should apply quite 
generally for any continuous mechanical instability of a nano- 
electromechanical system. 

Our result apply most directly to quantum dots situated 
on nanobeams or carbon nanotubes. Applying strain to the 
nanobeam in a controlled manner could, in principle, be ex- 
perimentally performed with the help of a break junction. In 
fact, it is quite conceivable that, e.g., some carbon nanotube 
structures happen to be close to the Euler instability due to 
specifics in the fabrication of individual nanostructures. Our 
predictions may be helpful to identify such "anomalous" (and 
potentially interesting) samples. 
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FIG. 9. Coordinate system used to describe the elastic properties of 
the rod. 



fundamental bending mode (n — 1) thus vanishes at the crit- 
ical force F c - Kq 2 , while all higher modes have a finite fre- 
quency and hence are neglected in what follows. For F > F c , 
the fundamental mode is unstable, and quartic corrections to 
the Lagrangian are necessary to ensure global stability. De- 
noting <jj\ — <d and X\ = X, expanding the Lagrangian (Al) 
to quartic order in the displacement and inserting the solution 
hi of the harmonic problem, we thus obtain the effective La- 
grangian close to the Euler instability, 



Appendix A: Elasticity theory of the Euler instability 

The elastic Lagrangian of a homogeneous rod of constant 
length L fixed at its two ends consists of three parts, 



L = -x 2 - 



L=T-<V h -<V F . 
The kinetic term reads 



T 



2 Jo 



dsh 2 , 



(Al) 



(A2) 



where cr is the linear mass density, s the arc length along the 
rod, and h(s, t) the displacement of the rod with respect to 
the w-axis (see Fig. 9). The bending energy, controlled by the 
bending rigidity k, is given by 



2 Jo 



ds 



df 


2 „ 


ds 


= 2 Jo 



ds 

1 -h' 2 



(A3) 



where f = («', h') is the tangent vector of the rod and primes 
denote derivative with respect to s. The last term in Eq. (Al) 
corresponds to the work done by the compression force F on 
the rod and reads 

<V F = -F(L - « max ) = -F f ds (l - Vl-AV 2 ) , (A4) 

where M max is the total extent of the rod along the M-axis (see 
Fig. 9). 

For small deflections (\h'\ <K 1), the Lagrangian (Al) be- 
comes, in harmonic approximation, 



-h' 



+ 2" 



1 



jo ^ 2 
with the corresponding Euler-Lagrange equation 

crh + ah"" + Fh" = 0. 



(A5) 



(A6) 



Equation (A6) can be solved by the eigenfunctions h(s, t) = 
Yj n h n {s, t) = Yin X„(t)g n (s), where g„(s) are the normal modes 
which follow from the solution of the characteristic equation. 
The frequency of the mode n reads 



2 K l( 2 F \ 

a> n =-q n [q n ~-) 



(A7) 



with q„ the associated wavenumber which depends on the con- 
sidered boundary conditions. The vibrational frequency of the 



with the effective mass 



f 

Jo 



-X L - -X\ 
4 



dsgf, 



and 



= J ds^'/Vi 2 - y^'i 4 )' 



(A8) 



(A9) 



(A10) 



with g\{L/2) = 1 such that X corresponds in Eq. (A8) to the 
actual displacement of the center of the rod. Notice that a pri- 
ori, a depends on the force F. However, close to the buckling 
instability, we can approximate F - F c in Eq. (A10). 

The parameters entering the effective Lagrangian (A8) and 
the corresponding vibrational Hamiltonian (2) are given in Ta- 
ble I for two types of boundary conditions: hinged end points 
(h\ 0L = h"\ 0L - 0) and clamped end points (h\ QL = h'\ = 
0). Notice that in the latter case, only an approximate solution 
of the Euler-Lagrange equation (A6) can be found, which is 
valid in the vicinity of the Euler instability, i.e., for F F c . 



Appendix B: Langevin dynamics of the mechanical degree of 
freedom 



For the convenience of the reader, we present a derivation 
of the Fokker-Planck equation (9). Our derivation is quite 
general as long as electronic transport is described by rate 
equations (sequential tunneling). We note that Fokker-Planck 
equations for nano-electromechanical systems appeared pre- 
viously, e.g., in Ref. 14. 

We assume that the electromechanical coupling takes the 
general form H c = /z(X)n(j, where h is an arbitrary function of 
the mode amplitude X. In the classical limit (fi\<A k%T), 
and in the sequential tunneling regime (KY <zz k^T), one can 
write a Boltzmann equation for the joint probability distribu- 
tion VniX, P, t) that the resonator is in charge state n (= 0, 1) 
and phase space point {X, P) at time t (Refs. 13, 44), 

d,P n = {HM - (-l)T i(X)Po + (-lyTioCXFi. (Bl) 

The Poisson bracket {/, g) = dxfdpg - dpfdxg describes 
the classical dynamics of the mechanical degree of freedom 
on the adiabatic potentials corresponding to the neutral and 
singly-charged states of the quantum dot with Hamiltonian 
r H„ = // v ib + h(X)n, where H v n, is the vibrational Hamilto- 
nian of Eq. (2). In Eq. (Bl), the position-dependent rates for 
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TABLE I. Parameters entering the effective Lagrangian (A8) and the vibrational Hamiltonian (2). 
Boundary conditions gi(s) F c or m 



»(") M'-B t r Ak1 

■*(") ■(?)' ;(t)'K) ¥ Mar 



tunneling of electrons in and out of the dot, roi(X) and Fio(X), 
respectively, account for the electronic dynamics. Notice that 
the following derivation does not depend on the specific form 
of these rates, and hence on the particular model for the quan- 
tum dot which one considers (e.g., metallic or molecular). 

Close to the Euler instability, the vibrational mode becomes 
slow ("critical slowing down") and the Poisson bracket in 
Eq. (Bl) can be considered as a small perturbation. If we ne- 
glect the Poisson bracket entirely in a first step, the stationary 
solution of Eq. (Bl) reads P = T 10 P/T, Pi = T 01 P/T, with 
P = P + Pi and r = F i + r 10 . Next, we account for the 
Poisson bracket perturbatively to leading order by making the 
ansatz 

P (X, P, t) = Y ^-P(X, P, t) - 5P{X, P, f), (B2a) 

1 (A) 

Pi (X, P, t) = ^-^-PiX, P, t) + 6P(X, P, t). (B2b) 
T(X) 

In the adiabatic limit, in which electronic tunneling is much 
faster than the vibrational dynamics (|o»| <K F), one would 
then expect that 6P is small compared to P itself. Inserting 
Eq. (B2) into Eq. (Bl), we obtain 



d,p = — d x p - F eS (X)d P p + (d x h) d P 6P 

m 



(B3) 



with the effective force F e fi (X) = -dxHmb ~ {dxh) rio(X) and 
no(X) = Vqi(X) /F(X) the average occupation of the dot for 
fixed position X, and 



d,5P = {<H , SP) - Y5P + {d x h) I -^^-d P P - 
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r 2 



dpSP 



P Eio^A-roi - r 01 d x r 



10 



p. 



(B4) 



m r 2 

So far, Eqs. (B3) and (B4) are exact, and we now make use 
of the adiabatic limit 6P <s P to solve them. In this limit, 
all terms in Eq. (B4) containing 5P are negligible compared 
to those involving P, except for the second term on the right- 
hand side which is multiplied by the total tunneling rate T » 
\oj\. We thus obtain from Eq. (B4) 

SP a (dxh) ^^d P P - - r ' oWl ~ r °'^ F ' V . (B5) 

Inserting this expression into Eq. (B3), we obtain the Fokker- 
Planck equation (9) (except for the purely extrinsic dissi- 
pative and diffusive parts proportional to the damping con- 
stant i] e , which can readily be obtained by coupling the me- 
chanical degree of freedom to a phonon bath), with damping 
coefficient rj(X) = - (dxh) (<9x«o) /T, and diffusion constant 
Z)(X) = 2(3 x /i) 2 »o(l-»o)/r. 



Appendix C: Current blockade for a harmonic oscillator 
linearly coupled to a SET 

In this Appendix, we detail the derivation of the Coulomb 
diamond and the resulting current blockade at low bias volt- 
age within the zero-temperature mean-field approximation of 
Sec. IV, in the specific case where the resonator coupled to 
the SET is purely harmonic. We follow and adapt Ref. 7 to 
the case of a single-level quantum dot, where the average oc- 
cupation of the island is given by Eq. (22). Specifically, we 
write the effective potential as 



with 



Veg(x) = v(x) + I dx'n (x'), 



u(x) — — - — (x — x) . 



(CI) 



(C2) 



Introducing x — x-x and t) g = v g — x, and using Eq. (22), the 
condition (21) for dynamical equilibrium at zero temperature 
reads 



• V"(X)X : 



1, X < Va — — , 

1 v u 

2' 5 H-2 <i<5 s + 2' 



(C3) 



0, 



v 



X > v„ + 



In order to solve Eq. (C3), one has to distinguish three cases 
(see Fig. 10): 

(i) v g > v/2: In that case, there only exists one solution to 
Eq. (C3) which corresponds to «o = 1 and thus is not 
conducting, jfi = -l/v"(x). Therefore, for v g > v/2, the 
system is always non conducting. 

(ii) |5g] < v/2: Here, two solutions can coexist: x\ (non 
conducting), and X\j% = —l/2v"(x) which corresponds 
to no = 1/2 and hence represents the conducting state 
of the effective potential. In order for the system to be 
conducting within our mean-field approximation, this 
latter solution must be the most stable. This happens 
whenever v > v+ = 2v g + 3/2v"(x). 

(iii) v g < -v/2: In this last case, an additional solution to 
Eq. (C3) exists on top of X\ and Jci/2, namely Xq — 
which is not conducting and corresponds to «o = 0. The 
conducting solution X1/2 1S the most stable one if v > v + 
and v > y_ = -2v g - l/2v"(x). 
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(iii) 



A„ \ 



(ii) 



(i) 



2v"(x) 



-l/2v"(x) 

FIG. 10. Sketch of the Coulomb diamond for a SET linearly coupled 
to a harmonic oscillator, delimited by the solid lines at v + = 2v % + 
3/2v"(x) and u_ = -2v g -l/2v"(x). The gray area indicates the region 
where current can flow. The dashed lines indicate the location of 
the Coulomb diamond without electromechanical coupling (F e = 0), 
delimited by v = ±2v g , which defines the three regions of the v-v g 
plane discussed in the text. 



The apex (defined by v+ = u_) of the resulting Coulomb dia- 
mond sketched in Fig. 10 is thus located at a bias voltage 



A„ = 



1 



2u"(x) 



(C4) 



which defines the energy gap below which current cannot flow 
through the system. The gate voltage corresponding to such a 
gap is given by v g = -l/2v"(x) + x. 



Appendix D: Analytical treatment of thermal fluctuations 

We present here a detailed analytical treatment of the be- 
havior of the current when one only considers thermal fluctu- 
ations in the Fokker- Planck equation (15), and compare it to 
the numerical results presented in Fig. 5 of Sec. V A. 

For temperatures much smaller than the gap, the onset of 
the current with bias voltage has a Fermi-function-like behav- 
ior (see dashed and dashed-dotted lines in Fig. 5). This be- 
havior can readily be checked by approximating the effective 
potential (25) for gate voltages at the apex of the Coulomb 
diamond [see Eq. (24)] by its zero-temperature expression, 
and by expanding the Boltzmann distribution (19) close to 
the corresponding minima of that potential. In this limit, the 
Boltzmann distribution is a superposition of weighted Gaus- 
sian peaks centered at these minima (except for |i5| <K JJa and 
x ^ 0, where the effective potential is purely quartic in x). For 
T <K A„, we find from Eq. (14) for the current 



In the vicinity of the instability (\6\ <K -^a), we have 



g(z) = 



1/3 2 2 / 3 F(5/2) /6A t ,\ 1/4 
2 2 / 3 r(5/2) /6A^ 1/4 

^ \ f 



o, 



<z < zi, 

Zl < Z < Z2, 

Z2 < Z < Z3, 
Z > Z3, 



(D3) 



where F(v) denotes the gamma function, and zi = (2 5 ^ 3 /3 - 
1)/(1 - 2-'/ 3 ), z 2 = 1/3(1 - 2- 1 / 3 ), and z 3 = 1/(1 - 2- 1 / 3 ). 
Notice that the discontinuities in Eqs. (D2) and (D3) are due 
to the fact that (meta)stable conducting or blocked minima of 
the effective potential are appearing or disappearing as one in- 
creases the bias voltage. Our approximate result (Dl) shows 
that at low temperature, the current below the gap is exponen- 
tially suppressed as a function of bias voltage. 

At larger temperatures, the current starts to be linear in the 
bias voltage, as can be seen from the dotted line in Fig. 5 for 
J/ A„ = 1. Indeed, expanding for T » \v g + v/2\ the current for 
fixed x, Eq. (13), as well as the Boltzmann distribution (19), 
we find using Eq. (14) 



I^C 



.eT v 
4 f 



(D4) 



with 



C 



I 



[ST 



at 3 



di// F (2/)exp|ytr - —y 



I 



-1 1st 2 

dyf F (-?/)expl —y 



af 3 



(D5) 



to first non-vanishing order in (u, \u g \)/T. For T » A„, we find 
C - 1 1 A, such that 



/ eT- 



eV 
l6k B T' 



(D6) 



which corresponds to the usual high-temperature current in 
absence of electromechanical coupling. We have checked that 
this result is in very good agreement with our numerical calcu- 
lation of Sec. V A. Of course, when the bias voltage becomes 
significantly larger than temperature, the current saturates to 
its maximal value I = eF/A. 



eT 
~A 



exp 



AT 



+ 1 



(Dl) 



in good agreement with our numerical results presented in 
Sec. V A. Far from the Euler instability (\6\ » y/a), the func- 
tion g(z) is given by 



g(z) 



2, < z < 2, 
0, z>2. 



(D2) 



Appendix E: Stationary solution of the Fokker-Planck equation 

We show here that in the adiabatic limit (ojo/F « 1) and 
for weak extrinsic dissipation (y e <s; 1), the stationary solu- 
tion of the Fokker-Planck equation (15) only depends on the 



ratio 



w /r' 



This behavior is examplified in Fig. 6 where the 
averaged current flowing through the nanobeam only depends 
on the latter ratio. 
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Introducing the new variables 

P 2 

E(x,p) = — + v eS (x), 

r dx 1 r x 
9{x > p)= \ \ 

Jx a A J Mi 



dx' 



y/2[E(x,p) - v eS (x')] 



(El) 
(E2) 



with E the energy of the mechanical degree of freedom and 6 
the time along a trajectory in phase space at a given E (xq is 
the initial position of the system for that particular trajectory), 
the Fokker-Planck equation (15) reads 



with dPst <s P st in Eq. (E3). To first order in (co /T, y e ) <sc 1, 
we obtain 



0=- — +Wj ) + re )^ + p 



d(x) fMdPgt 2 d 2 P st 



(E5) 



Averaging this equation over one period T of the motion in 
phase space, and using the periodicity of 6P st in 8, we have 



dP 



dP , , s / 2 dP d8dP\ 

(d(x) ~\(dP 2 d 2 P 8 2 8dP 
+ [— +7 ' T j[dE +P dF 2+ di 2 ^ 



d8\ 2 d 2 P 88 d 2 P \ 
dp)~d8~ 2+ P d~pdEd8j 



(E3) 



Noticing that the current-induced fluctuation and dissipation 
are both proportional to coq/T [cf. Eqs. (17) and (18)], we ob- 
tain that for wo/r = y e — 0, the stationary solution of Eq. (E3), 
d T Pst = 0, is independent of 8. This suggests to insert the 
ansatz 



P st (E,G) =P st (E) + 6P st (E,8) 



(E4) 



d(x) 



d \\ dP s , 



where 



(f{x,p)) E =7f ^ d6f(x,p). 



■T 



(El) 



Since y(x) and d(x) only depends on coq/T via their prefac- 
tors, it is clear from Eq. (E6) that the stationary solution of 
the Fokker-Planck equation only depends on a>o, T, and y e 
through the ratio . 



1 R. G. Knobel and A. N. Cleland, Appl. Phys. Lett. 81, 2258 

(2002) . 

2 R. G. Knobel and A. N. Cleland, Nature (London) 424, 291 

(2003) . 

3 A. D. Armour, M. P. Blencowe, and K. C. Schwab, Phys. Rev. 
Lett. 88, 148301 (2002). 

4 L. Y. Gorelik, A. Isacsson, M. V. Voinova, B. Kasemo, R. I. 
Shekhter, and M. Jonson, Phys. Rev. Lett. 80, 4526 (1998). 

5 F. Pistolesi and R. Fazio, Phys. Rev. Lett. 94, 036806 (2005). 

6 O. Usmani, Ya. M. Blanter, and Yu. V. Nazarov, Phys. Rev. B 75, 
195312 (2007). 

7 F. Pistolesi and S. Labarthe, Phys. Rev. B 76, 165317 (2007). 

8 J. Koch and F. von Oppen, Phys. Rev. Lett. 94, 206804 (2005). 

9 J. Koch, F. von Oppen, and A. V. Andreev, Phys. Rev. B 74, 
205438 (2006). 

10 R. Leturcq, C. Stampfer, K. Inderbitzin, L. Durrer, C. Hierold, E. 
Mariani, M. G. Schultz, F von Oppen, and K. Ensslin, Nat. Phys. 
5, 327 (2009). 

11 D. Mozyrsky, M. B. Hastings, and I. Martin, Phys. Rev. B 73, 
035104 (2006). 

12 F. Pistolesi, Ya. M. Blanter, and I. Martin, Phys. Rev. B 78, 
085127 (2008). 

13 A. D. Armour, M. P. Blencowe, and Y. Zhang, Phys. Rev. B 69, 
125313 (2004). 

14 Ya. M. Blanter, O. Usmani, and Yu. V. Nazarov, Phys. Rev. Lett. 
93, 136802 (2004); 94, 049904(E) (2005). 

15 C. B. Doiron, W. Belzig, and C. Bruder, Phys. Rev. B 74, 205336 
(2006). 

16 G. A. Steele, A. K. Huttel, B. Witkamp, M. Poot, H. B. Meer- 
waldt, L. P. Kouwenhoven, and H. S. J. van der Zant, Science 



325, 1 103 (2009). 

17 B. Lassagne, Y. Tarakanov, J. Kinaret, D. Garcia-Sanchez, and A. 
Bachtold, Science 325, 1107 (2009). 

18 L. Euler, in Leonhard Eider's Elastic Curves, translated and anno- 
tated by W. A. Oldfather, C. A. Ellis, and D. M. Brown, reprinted 
from ISIS, No. 58 XX(1), 1744 (Saint Catherine Press, Bruges). 

19 L. D. Landau and E. M. Lifshitz, Theory of Elasticity (Pergamon 
Press, Oxford, 1970). 

20 M. R. Falvo, G. J. Harris, R. M. Taylor II, V. Chi, F. P. Brooks Ir, 
S. Washburn, and R. Superfine, Nature (London) 389, 582 (1997). 

21 S. M. Carr and M. N. Wybourne, Appl. Phys. Lett. 82, 709 (2003). 

22 S. M. Carr, W. E. Lawrence, and M. N. Wybourne, Europhys. Lett. 
69, 952 (2005). 

23 D. Roodenburg, J. W. Spronck, H. S. J. van der Zant, and W. J. 
Venstra, Appl. Phys. Lett. 94, 183501 (2009). 

24 S. M. Carr, W. E. Lawrence, and M. N. Wybourne, Phys. Rev. B 
64, 220101(R) (2001). 

25 P. Werner and W. Zwerger, Europhys. Lett. 65, 158 (2004). 

26 V. Peano and M. Thorwart, New J. Phys. 8, 21 (2006). 

27 S. Savel'ev, X. Hu, and F. Nori, New J. Phys. 8, 105 (2006). 

28 G. Weick, F. Pistolesi, E. Mariani, and F. von Oppen, Phys. Rev. 
B 81, 121409(R) (2010). 

29 A. K. Huttel, G. A. Steele, B. Witkamp, M. Poot, L. P. Kouwen- 
hoven, and H. S. J. van der Zant, Nano Lett. 9, 2547 (2009). 

30 Typical parameters for carbon nanotubes of length L and 
radius R are: k = 6.3i? 3 [nm]/jeV.m, cr = 4.2 X 
10" 15 i?[nm]kg.m _1 , such that m = 1.6 x 10 2[ R[nm]E\jjm] kg, 
io = 610fl[nm]Zr 2 [jum]MHz, F c = 40fi 3 [nm]Zr 2 [jum] pN, and 
a = 1.5 X 10 27 # 3 [nm]Zr 5 [pm] eV.nT 4 . 

31 For simplicity we consider here spinless fermions. Within the ap- 



15 



proximations of this paper, the spin would add only a multiplicity 
factor in the tunneling rates. 

32 E. Mariani and F. von Oppen, Phys. Rev. B 80, 15541 1 (2009). 

33 S. Sapmaz, Ya. M. Blanter, L. Gurevich, and H. S. J. van der Zant, 
Phys. Rev. B 67, 235414 (2003). 

34 W. Izumida and M. Grifoni, New J. Phys. 7, 244 (2005). 

35 S. Sapmaz, P. Jarillo-Herrero, Ya. M. Blanter, C. Dekker, and H. 
S. J. van der Zant, Phys. Rev. Lett. 96, 026801 (2006). 

36 K. Flensberg, New J. Phys. 8, 5 (2006). 

37 More precisely, if gX 2 n A /2 is the intrinsic coupling term (see 
Ref. 28), we can neglect this term with respect to the linear one if 

g « F 2/ V' 3 . 

38 R. Hussein, A. Metelmann, P. Zedler, and T. Brandes, Phys. Rev. 
B 82, 165406 (2010). 

39 U. Weiss, Quantum Dissipative Systems (World Scientific, Singa- 
pore, 1993) 

40 P. Mohanty, D. A. Harrington, K. L. Ekinci, Y. T. Yang, M. J. 
Murphy, and M. L. Roukes, Phys. Rev. B 66, 085416 (2002). 

41 C. Seoanez, F. Guinea, and A. H. Castro Neto, EPL 78, 60002 
(2007). 

42 C. Seoanez, F. Guinea, and A. H. Castro Neto, Phys. Rev. B 77, 



125107 (2008). 

43 C. Seoanez, F. Guinea, and A. H. Castro Neto, Phys. Rev. B 76, 
125427 (2007). 

44 F. Elste, G. Weick, C. Timm, and F. von Oppen, Appl. Phys. A 93, 
345 (2008). 

45 Notice that at equilibrium (V = 0), the fluctuating and dissipa- 
tive parts of the Fokker-Planck equation (9) fulfills the fluctuation- 
dissipation theorem, i.e., [D(X)/2 + ri s k B T]/[n(X) + n c ] = k B T. 

46 Notice that in actual experiments on suspended carbon nanotube 
quantum dots, the ratio ti)<j/r is typically very small, of the order 
of 10~ 2 (seeRefs. 16, 17, 35). 

47 G. Weick et ah, unpublished. 

48 P. M. Chaikin and T. C. Lubensky, Principles of Condensed Mat- 
ter Physics (Cambridge University Press, Cambridge, 1995). 

49 At equilibrium (v = 0), the effective temperature (26) corresponds 
to the temperature of the electron reservoirs T, cf. Ref. 45. 

50 Notice, however, that this does not imply a heating of the elec- 
tron reservoirs as these are considered to be macroscopic. Thus, 
there is no back-action of the mechanical degree of freedom on 
the electronic leads which remain at the equilibrium temperature 
T. 



